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We explore the possibility of inducing a topological insulator phase in a honeycomb lattice lack- 
ing spin-orbit interaction using a metallic (or Fermi gas) environment. The lattice and the metallic 
environment interact through a density-density interaction without particle tunneling, and integrat- 
ing out the metallic environment produces a honeycomb sheet with in-plane oscillating long-ranged 
interactions. We find the ground state of the interacting system in a variational mean-field method 
and show that the Fermi wave vector, kp, of the metal determines which phase occurs in the honey- 
comb lattice sheet. This is analogous to the Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism 
in which the metal's &f determines the interaction profile as a function of the distance. Tuning 
fci? and the interaction strength may lead to a variety of ordered phases, including a topological 
insulator and anomalous quantum-hall states with complex next-nearest-neighbor hopping, as in 
the Haldane and the Kane-Mele model. We estimate the required range of parameters needed for 
the topological state and find that the Fermi vector of the metallic gate should be of the order of 
37r/8a (with a being the graphene lattice constant). The net coupling between the layers, which 
includes screening in the metal, should be of the order of the honeycomb lattice bandwidth. This 
configuration should be most easily realized in a cold-atoms setting with two interacting Fermionic 
species. 



I. INTRODUCTION 

In 2005 Kane and Mele^l proposed that graphene, a 
one atom thick graphite can exhibit helical edge states 
which arc protected against weak perturbations by topol- 
ogy. This idea stems from the Haldane modeP of a quan- 
tum Hall state without magnetic field. This special topo- 
logical phase would have emerged as the ground state 
of graphene if the intrinsic spin-orbit coupling was large 
enough and in particular, larger than the Rashba spin- 
orbit coupling. Unfortunately, ab-initio calculations^ 
have found that this requirement is far from being ful- 
filled in graphene and the desired intrinsic spin-orbit cou- 
pling is smaller than O.OOlmeV (about 0.01K). Neverthe- 
less, this direction has led the way to a wealth of theoreti- 
cal predictions^*! of topological states. In parallel, exper- 
iments have show n that robust helical edge states exist 
in t hree d imension ;] 10 ! 11 ! and in two dimensional quantum 
wellsPEal 

Despite the fact that Graphene's intrinsic spin-orbit 
coupling may be too weak to produce the desired topo- 
logical insulator phase the hope to achieve such a phase 
in the honeycomb lattice has not died. In this work 
we explore the possibility of engineering a honeycomb- 
lattice-based or Graphene-based topological insulator by 
employing an appropriate environment. Our results are 
relevant either to a Graphene with a metallic gate, or to 
a Fermi-gas residing in a honeycomb optical lattice (see, 
e.g., Refp31), and interacting with a second Fermi-gas of 
similar (two-dimensional) density that does not couple to 
the optical lattice. Note that another path for producing 
topological states in optical lattices was explored in Ref. 
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An important ingredient for topological behavior is a 



bulk gap that arises through avoided level crossing be- 
tween two bands and a resulting band inversion. A half 
filled honeycomb lattice (e.g., undoped graphene) has 
band crossing at the two Dirac points^^ but without 
any band gap. An opening of a gap in the bulk is possible 
in a few ways, for example by substrate effects^ or a large 
lattice distortion!^. However, not all gaps are alike. In a 
two valley model like graphene, non-trivial topology can 
only occur if the gap function changes sign between the 
two valleys. This sign change can occur, for instance, due 
to intrinsic spin-orbit coupling. On the lattice, this cou- 
pling may be approximate by imaginary hopping between 
next nearest neighbor sites and this yields the Kane-Mele 
model. The sign of the imaginary hopping depends not 
only on the direction of the hopping but also on whether 
the path from one site to another contains a left or right 
turn. In momentum space, the imaginary hopping takes 
the form of a function that changes sign between Dirac 
valleys, and between spin polarizations. 

Apart from intrinsic spin-orbit coupling, this type of 
directional hopping term has been found to aris e from an 
interaction between next nearest neig hbor sitesPESl i t 
has also been found that interactions may lead to a topo- 
logical insulator in a decorated honeycomb and Kagome 
lattices^. 

The main obstacle for realizing an interaction-driven 
topological insulator is the required interaction profile in 
real space. In a honeycomb lattice, for example, this 
amounts to an interaction which is strongest on next- 
nearest-neighbor bonds. When this requirement is not 
fulfilled, other phases (density waves, lattice distortions, 
etc.) may occur. An on-site repulsion U (arising due to 
Coulomb interactions in electrons, or on-site scattering 
for cold atom), enhances the tendency for magnetic or- 
dering. The next strongest interaction (for electrons) is 
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the nearest- neighbor repulsion, V\, which tends to favor 
charge/spin density waves. A second-nearest-neighbor 
interaction V2, however, enhances precisely the propen- 
sity to induce an imaginary hopping term leading to the 
topological phase. It is hard to find, however, lattice 
systems where the second-nearest neighbor interaction 
trumps the shorter term repulsions; therefore, the topo- 
logical phase is not likely to arise due to intrinsic inter- 
actions. In addition, if V\ and V2 are of similar order, 
a Kekule bond order may occur^. Extensive studies on 
interacting graphene have shown a variety of competing 
phaseJ22HSSl. Below we discuss the effects of producing 
the longer range interactions required for topological or- 
der using the Friedel oscillations and RKKY interactions 
occurring in a metallic (Fermi-gas) environment. Our 
analysis considers the possibility of forming both topo- 
logical phases as well as other competing phases. 



II. TUNING THE INTERACTION 

In this work we propose a setup in which a 'knob' can 
be turned to control the relative ratio between the dif- 
ferent interaction terms. Our idea is motivated by the 
RKK ^25H!3 interaction in metals. The original RKKY 
model was written to describe how impurity spins in- 
teract with each other over a long distance through a 
medium of electronic states. The interaction does not 
simply fall as a power law of the distance between the 
spins like the usual Coulomb interaction but instead os- 
cillates. The oscillations can be viewed as Friedel oscilla- 
tions in the metal caused by the impurity spins. The in- 
teraction profile depends on the polarization of the metal 
and its functional form is — cos(2kpr)/r 3 in three dimen- 
sions and — cos(2kpr)/r 2 in two. The Fermi gas back- 
ground is considered as a simple parabolic band. 

With this scenario in mind we propose to generate in- 
teractions in the honeycomb lattice that are mediated 
by an additional Fermi gas or metallic gate. An elec- 
tronic realization of such a system consists of a metallic 
layer which is put above a honeycomb lattice such that 
the quasiparticles in it interact with the metallic-gate 
quasiparticles without direct hopping between the two 
systems. A more promising avenue, which our analysis 
more readily addresses, is a cold-atoms realization, con- 
sisting of two different species of Fermions, only one of 
which is trapped in an optical honeycomb lattice, with a 
contact interaction between the two species. The effect 
of the coupling between the honeycomb fermions and the 
free gas can be found by integrating out the metallic de- 
grees of freedom. This results in the honeycomb fermions 
acquiring a long range interaction whose magnitude de- 
pends on the polarization of the free Fermi gas. For sim- 
plicity, we will assume that the fermions are electron-like, 
and refer to them as electrons and fermions interchange- 
ably below. 



III. FORMALISM 

Our starting point is electrons hopping on a honey- 
comb lattice and interacting through an on-site density- 
density interaction with fermions in a metallic layer: 
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where c ia creates an electron with spin a on the honey- 
comb lattice site i, d^ creates an electron with spin a in a 
Bloch state k of the metal and n\ and nf are the number 
operators at the position i of the honeycomb lattice and 
metal electrons respectively. In the metal this should be 
understood as a coarse grained density integrated over a 
small area, A e tt around the position i. This area should 
be estimated by the screening in the metal, or in the cold- 
atom case, by the probability density of the two species. 
We set it to the lattice constant a in order to simplify 
the notation. A larger area may enhance the interaction 
but will also make the approximate interaction term in 
Eq. [I] less accurate. We estimate the value of a 2 in the 
discussion section and in the appendix. 

In the path integral language the action is given by 
a time integration over the Lagrangian which is derived 
from the above Hamiltonian and the integration over the 
metallic Fermion operators is possible due to their 
quadratic form. The result is an interacting theory of 
quasiparticles on the graphene sheet with their action 
given by: 



dr 



(ij) a 
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where the integration is over the imaginary time variables 
r and t' and Tlij(r — r') is the polarization operator of 
the metal. The interaction above is long ranged and dy- 
namic. Its functional form in momentum and Matsubara 
frequency space is given by the Matsubara sum of the 
polarization bubble of the Fermi sea: 



n(ifi,q) 
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where np is the Fermi-Dirac occupation number, which 
includes the chemical potential. This polarization opera- 
tor leads to a long ranged, time dependent interaction in 
the honeycomb lattice. In the current work, for simplic- 
ity, we explore the limit of static interaction. The static 
limit is valid if the velocity in which information travels 
in the honeycomb lattice is much slower than the typi- 
cal velocity in the metal. The comparison can be made 
in terms of energy. In the graphene sheet we take the 
energy of the gap that opens as a result of the topolog- 
ical order and in the metal it is the Fermi energy. We 
found that there is indeed a range of parameters where 
the developed gap is small compared to the metal's Fermi 
energy. 



(a) U = 



(b) U = t 



(c) U = 2t 



FIG. 1: Phase diagram for graphene with induced long ranged interaction of the form of Eq. |4| The x-axis is 2kp of 
the metallic layer (measured in units of l/a where a — 2A6A is the graphene lattice constant) which governs the 
relative strength of interaction in different distances and the y-axis is the overall interaction strength which is 
determined by the coupling a. The on-site interaction parameter U was scanned. Here we show (a) U = 0, (b) 
U = t and (c) U = 2t where t is the hopping amplitude. The different symbols (colors online) represent different 
phases with: circle (black) semi-metal, square (blue) charge density wave, diamond (gray) spin density wave, 
up-triangle (purple) superconductor, down-triangle (orange) Kckule, open circle (red) anomalous Hall, open square 

(green) anomalous spin Hall. 



IV. METHOD AND RESULTS 

The static limit of the polarization bubb le is the fa- 
mous Lindhard function in two dimension d 28 * 29 !. It is 
given by: 

n(»^o,Q = - fl2 cos(2 9 fcFr) (4) 

where r is the distance between the interacting sites, &f 
is the Fermi wavevector of the metal, the a 2 factor rep- 
resents the effective interaction area A e f f and 1 jt is the 
estimated metallic density of states at the Fermi level. 
In order to make a crude estimate of the effective inter- 
action strength in the layer we assume that the interac- 
tion between the electrons on the graphene sheet and the 
metal electrons is e 2 /d where the distance d is of the or- 
der of 10 A. This gives a coupling a which is of the order 
of the bandwidth t (2-3eV). However, when screening is 
taken into account (see the appendix) the bare coupling 
is reduced by a factor of exp(— 4k pd) and therefore the 
required distance is much smaller. This poses a great 
challenge since at small distances tunneling may occur. 
We present results (Fig. [lj for this order of interactions 
where (a/t) 2 is scanned from to 4. 

The standard mean field decomposition leads to self- 
consistency equations that are usually solved by itera- 
tions. The interaction terms are decomposed in the vari- 
ous channels. For example, the next-nearest neighbor in- 
teraction has been decoupled by Raghu et al. into the de- 
sired second neighbor hopping interaction: V^Cicjcj — > 

y2(Xij c j c i + h.c—\xij\ 2 ) and the self-consistency equation 

reads Xij = ( c l c j)- This decomposition is reasonable, 
however, one should take care to decouple all interaction 
terms in all channels. When this is not done correctly 
the combined effect of multiple terms may result in errors 
and phases may be missed. Given the drawbacks of the 
standard mean-field decoupling we adopt the variational 
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FIG. 2: An example of order parameter magnitudes as 
a function of the metal's Fermi wavelength (measured 
in units of l/a where a = 2AA is the graphene lattice 
constant). The data is given for a constant U = and 
Vq = 4t. The curves (color online) represent the order 
parameter magnitude in units of the bare hopping t 
with: (blue) square - charge density wave, (red) circle - 
anomalous Hall and (purple) triangle - 
superconductivity. 



mean field approach^. This is a systematic procedure 
which does not suffer from the above problems. We gen- 
erate a variational wavefunction that is a solution to a 
quadratic auxiliary Hamiltonian. This generating Hamil- 
tonian, Hg en , contains the graphene nearest neighbor 
hopping and a variety of order parameters. These include 
charge and spin density waves in which the two sublat- 
tices have different charge/spin density; the Kekule dis- 
tortion, which is a structural distortion known to occur in 
carbon nanotubes^ and may have interesting topological 
excitations^!, the anomalous Hall or spin Hall state and 
superconductivity. Superconductivity is included here 
since at certain Fermi vectors some interaction coeffi- 
cients may be attractive. The resulting superconducting 
state has an interesting momentum-dependent order pa- 
rameter which resembles the graphene's band dispersion. 
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Near the valleys the superconducting order parameter 
has a p x iiPy form with the sign changing between the 
two valleys 33 . Another interesting suggestion for this re- 
gion is Kekule superconductivity^. 

The generating Hamiltonian is a 12 x 12 matrix for each 
momentum k with the following structure. 



^gen(k) 



tPf g 
$ -hi 



hf = 



f Haa Hab 
\Hba Hbb 



9 = 



( Gaa Gab 
\Gba Gbb 



(5) 



where the 2x2 structure of H gen is the Nambu space 
with h containing the usual particle-hole terms and g 
containing the pairing amplitude. The matrices h and g 
are further split into the A and B sublattices and due 
to the Kekule distortion (which enlarges the unit cell to 
include three atoms of each sublattice) the matrices H nm 
and G nm of dimension 3. 



Haa = S] 



Hab 

Hba 
Gaa 

Gab 
Gba 




t 3 e 



ika 2 



c ika 3 \ 

ika 2 



AB 



(6) 



Here p^ and p^ are the order parameters for density 
waves in the up and down spin respectively, <Skt/i 



ik-ai 



Q ika 2 



s lk ' a3 ) are the Fourier transform 



of the second nearest neighbor hopping with A^ and A^ 
being the order parameter for topological order in the up 
and down spin and the vectors a! , a 2 and a 3 are (hexag- 
onal) lattice vectors. (a x = ax, a 2 / 3 = a(^x ± ^y)). 
Following Weeks and FranJ^Hl we parametrize the lattice 
distortion by tj = t + St + rjj where St is a uniform shift 
and the Kekule texture is given by rjj — r)cos(^j + 4>) 
with r] being the order parameter and 4> a free parameter 
which does not affect the size of the gap. The down spin 
part of the generating Hamiltonian is obtained from 
h\ by replacing p^ — > and Af — ^ A4.. 

The trial wavefunction, which is the ground state of 
the Hamiltonian in Eq. [5] is used to calculate the full 
interacting variational energy. This is done by calculating 
the expectation of the Hamiltonian: 



" = -*£$: 



it n 4 + 2_\ V s nin i+5 (7) 
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where U represents an effective on-site interaction and V$ 
represents the interaction between two sites separated by 



the distance S which is taken up to the seventh neighbor 
and is given by the coupling a 2 times the RKKY inter- 
action in Eq. |4j The on-site interaction is estimated at 
3.5eV (about 1.3tp^but since the uncertainty is large, we 
choose to scan different values of U. The grid size we use 
is 64 x 64 unit cells (of six atoms each) . In order to save 
computing time we calculate the correlation functions in 
momentum space, Fourier transform them and use the 
resulting real-space correlations to evaluate the interac- 
tion energy for any given distance. The next step is to 
minimize the energy with respect to the order parame- 
ters and the minimum determines the mean-field ground 
state. We use the standard downhill simplex method 3 -^ to 
do this minimization. The results are shown in Figs. [l][2] 



V. SUMMARY AND DISCUSSION 

In this Paper we analyzed the effect of a metallic en- 
vironment on fcrmions residing in a honeycomb lattice. 
The main idea is that the environment, through its Lind- 
hard response, can induce long range interactions in the 
honeycomb lattice, which are conducive to the formation 
of topological phases, such as the anomalous and spin- 
Hall states. Our analysis consisted of a variational study 
of the many possible orders of a honeycomb lattice in the 
presence of an oscillating long-range interaction. To test 
our analysis, and to make a connection with previously 
derived results we used an interaction cut-off at the next 
nearest neighbor length and plot a phase diagram for a 
spinless V\ — V% model. Our results are similar to those 
obtained by Weeks and Franz^l when spinless Fermions 
are considered. It is interesting to note that the Kekule 
phase vanishes as soon as the spin and on-site interaction 
are introduced. 

When we introduce the oscillating long-range compo- 
nent of the interaction, we indeed find that topological 
phases, among other phases, may be induced in a hon- 
eycomb lattice through interaction with a metallic envi- 
ronment. In particular, topological order occurs in the 
case of strong coupling between the two species, when 
the Fermi wave number of the metal is between n/Aa 
and 7r/2a (where o is the lattice constant of the honey- 
comb). This is shown in Fig. [l] where the anomalous 
Hall and spin Hall effects are marked by open symbols. 

The most promising path for realizing our proposal is 
in a cold-atoms contexts. A cold-atoms realization would 
consist of using two species of Fermionic gases. The first 
would be confined to a honeycomb lattice (and could have 
one or two hyperfine states to imitate spinless or spin- 
half particles). The second species would be confined to 
the honeycomb lattice plane, but without being sensitive 
to the optical lattice that traps the first species. This 
setup allows the realization of contact interactions be- 
tween the two species through a Feshbach resonance, as 
well as the on-site intraspecies interactions described in 
Eq. without significant limitations. 

A realization of the scheme using electrons would em- 



5 



ploy a gated graphene layer. Unfortunately, it suf- 
fers from significant drawbacks, since it is a long- 
range Coulomb interaction which couples between the 
Graphene and the gate. First, when U = 0,t (see Fig.JT]- 
a,b), the coupling constant a 2 that is needed in order to 
induce topological order is about 2t 2 . Defining d as the 
distance between the layers, the interaction parameter 
is a = exp(— 2fcp<i). The exponential suppression 

is the result of convolving the density modulations in 
the gate with the Coulomb potential as discussed in Ap- 
pendix [a] Even choosing a gate as close as d = 6 A yields 
a ~ t/50. 

In addition, the large wavevector needed in the gate 
corresponds to a large density of the two dimensional 
electron gas which is very difficult to achieve in met- 
als or semiconductors. Perhaps this could be mitigated 
by using doped graphene as the metallic layer to in- 
duce the interactions. At low density, the Friedel os- 
cillations and resulting RKKY interactions in graphene 
are the result of scattering in the vicinity of the same 
valley point, i.e., short wavevector, and between differ- 
ent valleys. While intcr-vallcy scattering is suppressed 
due to the Klein paradox and decay as 1/r 3 , the inter- 
valley scattering is allowed and is of the right order of 



magnitude for our purposes' 37 * 38 *. This suggests that the 
ideal gate for the realization of topological phases (if suf- 
ficient proximity could be achieved) is another layer of 
graphene that is rotates by 30° with respect to the bot- 
tom one and is doped such that the inter-valley scatter- 
ing vector matches the ideal wave-vector. The challenge 
then becomes to prevent the tunneling between the two 
rather close Graphene planes. Note that a similar setup 
was considered in RefsP^ED 
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Appendix A: Estimate of the induced interaction in the graphene layer 



In this section we concentrate on the electronic realization of our proposal, and consider the details of gate-induced 
interaction between two points on a graphene sheet, r[ and r' 2 , whose distance is i?, taking into account screening, and 
the distance d to the gate. We use the known result for the interaction between these points and two points on the 
metallic layer r± and r 2 which are close to r[ and r 2 but not necessarily the same. The interaction term is therefore: 



e 4 



[w, /^^"'^""^y^y^) (ad 

J J ' \n - r 2 \\r 2 - r' 2 \ 



We would like to integrate out the metal by integrating over the positions r,: 

« e 4 / dV^rVKKfa) / d 2 nd 2 r 2 ^ m ^) nm ^)) 



In - r' 2 \\r 2 - r' 2 \ 
d 2 r{d 2 r' 2 n^r[)n^r' 2 )V eSf {r' i y 2 ) (A2) 

The effective interaction can be calculated using the polarization of the metal: 

1 fj2_ * (n m (n)n™(r 2 )) 



\t\ - r' x \\ri - r' 2 



4n,M 1 f ,2 ,2 COS(2fc F (pi - p 2 )) 1 
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= -e 4 D(0)4 I #Pi<Pp2 ~Z)"'"f yri — / - / ( A3 ) 

{ 'a 2 J n ^ 2 (2k F ( Pl - p 2 )) 2 ^d 2 + (p 1 ~p' 1 ) 2 ^/d 2 + ( P 2-p' 2 ) 2 ^ ' 

where D(0) = l/2ta 2 is the metallic density of states at the Fermi level and a ~ 2AA is the lattice constant. We 
have explicitly written the distances |ri — r[\ and \r 2 — r' 2 \ taking into account the distance between the layers d and 
denoting the distance in each layer by p. We rewrite the effective interaction: 

v «< - W^ ^^ylffrf'x* + a <A4) 

where we measure p\ from p\ and p 2 from p' 2 in the plane. With this choice x = 

yjR 2 + pi + p\ — 2R{pi cos(£?i) — p2 cos{9 2 )) — 2p\p 2 cos(6*i — 6%) is the distance between p\ and p2- Now, we 
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assume that R = \p[ — p' 2 \ is larger than px,p2 hi the relevant part of the integral (basically assuming \R\ <C d) and 
we can simplify the expression. We also replace x 1 in the denominator by R 2 and linearize in r/R in the cosine. This 
gives: 

e 4 f f cos(2fc F ( J R-piC0s(6»i) + p 2 cos(6 2 ))) 



2ta\2k F R) 2 J P1 P1 J W W v^+PiH^+Pl) 

J$(2k F pi)J (2k F p 2 ) cos(2k F R) 
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2ta±{2k F RY \j r ^(^2+^ 
2 e 4 cos(2kpR) 



2ta 4 (2k F R) 2 \ 2k F 
_ 2 e 4 cos(2fc F fl) 4fcFd 

_ 77 8ta 4 fc 4 i? 2 6 (A5j 
where we have used f dflcos^+frcos^)) = 27rJ (6) cos(a) in the angular integrals. Indeed, this expression is suppressed 

(2fc F a) 2 



(2n) 

by 7^ \ 2 exp(— 4k F d) with respect to the bare coupling square. 
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